Boltzmann theory of engineered anisotropic magnetoresistance in (Ga,Mn)As 
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We report on a theoretical study of dc transport coeffi- 
cients in (Ga,Mn)As diluted magnetic semiconductor ferro- 
magnets that accounts for quasiparticle scattering from ion- 
ized Mn 2+ acceptors with a local moment S — 5/2 and 
from non-magnetic compensating defects. In metallic samples 
Boltzmann transport theory with Golden rule scattering rates 
accounts for the principle trends of the measured difference 
between resistances for magnetizations parallel and perpen- 
dicular to the current. We predict that the sign and magni- 
tude of the anisotropic magnetoresistance can be changed by 
strain engineering or by altering chemical composition. 

In most (III,V) semiconductors, Mn 2+ substitution on 
a cation (column III element) site introduces an S = 5/2 
local moment and a valence band hole Q). Mn^IIIi-^V 
diluted magnetic semiconductors (|-|| (DMSs) are fer- 
romagnetic and metallic for Mn fractions larger than 
x ~ 1%. Many magnetic properties of the most ro- 
bustly ferromagnetic samples, those that have ~ 5% or 
more Mn and are annealed to reduce the density of com- 
pensating defects, appear to be adequately explained by 
virtual crystal approximation models in which disorder 
is ignored For example in bulk DMSs, this ap- 

proach can account for ferromagnetic critical tempera- 
tures ^100 K the correlation between magneto- 
crystalline anisotropy and substrate lattice constant , 
the size and sign of the anomalous Hall effect and 
several optical properties Q. It has also been used to 
describe properties of DMS heterostructures [p~2| p~4|| for 
which the simplifications afforded by neglect of disorder 
are particularly helpful. 

In this letter we investigate theoretically the T — dc- 
transport coefficients of (Ga,Mn)As ferromagnetic semi- 
conductors. We find that relaxation-time-approximation 
solutions of the Boltzmann equation provide anisotropic 
magnetoresistance (AMR) estimates that are in good 
agreement with experiments 16 . Our results suggests 
that transport properties of these metallic ferromagnets 
can be understood within a conventional framework in 
which disorder is treated as a weak perturbation. We find 
that the conductivity varies by several percent when the 
magnetic order parameter is reoriented by a weak mag- 
netic field, and predict that the magnitude and sense of 
this change depends on the chemical composition and on 
the substrate on which the thin film DMS ferromagnet is 
epitaxially grown. This spontaneous magnetoresistance 
anisotropy is the transport analog of magneto-crystalline 
anisotropy |t],|8| which has approximately the same size 



relative to the total condensation energy of the ordered 
state. All results presented in this paper are for the 
(Ga,Mn)As DMS's. A large database that details our 
predictions for the AMR of many other host semicon- 
ductors over a wide range of compositions and strains is 
available on the internet ]i7|j . 

We consider a microscopic Hamiltonian in which va- 
lence band holes interact with randomly located spins 
of substitutional Mn 2+ impurities via exchange interac- 
tions, and with randomly located ionized defects and each 
other via Coulomb interactions. Focusing on T — 0, we 
assume that the Mn spins are fully aligned in the ferro- 
magnetic ground state. In the virtual crystal approxima- 
tion, the interactions are replaced by their spatial aver- 
ages, so that the Coulomb interaction vanishes and hole 
quasiparticles interact with a spatially constant Zeeman 
field. The unperturbed Hamiltonian for the holes then 
reads Hq = Hl + J p dN Mn 2+ SCl ■ s , where Hl is the six- 
band Kohn-Luttinger Hamiltonian Q, O is the Mn local 
moment orientation, J p d = 55 meV nm 3 Q is the local- 
moment - valence-band-hole kinetic-exchange coupling 
constant, N Mn 2+ is the density of ordered Mn local mo- 
ments, and s is the envelope-function hole spin operator 
We use the relaxation-time-approximation solution 
to the semiclassical Boltzmann equation to estimate the 
dc conductivity tensor: 
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where T n g is the quasiparticle elastic scattering rate, n 

and k are the band and wavevector indices of the va- 
lence band Bloch states of the unperturbed system, and 
E g are the spin-split band energies of the ferromag- 
netic state. In Eq. (|l|) we have omitted the asymmetric 
terms in the off-diagonal elements of a a p that contribute 
to the anomalous Hall conductivity, discussed in detail 
elsewhere pxj] . The symmetric off-diagonal elements, de- 
scribed by Eq. (|l]), vanish when the magnetization is 
aligned along one of the cube edges of the host lattice. 

In our model, itinerant holes are scattered on substi- 
tutional Mn 2+ impurities by a Thomas-Fermi screened 
Coulomb potential and by a magnetic-kinetic-exchangc 
potential. For majority-spin holes both potentials are 
attractive while for minority-spin holes the magnetic po- 
tential becomes repulsive. We estimate the transport 
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weighted scattering rate from Mn 2+ impurities using 
Fermi's golden rule: 



T n,k ~ Y Mn2+ 2^ 
n' 
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where the scattering matrix element was approximated 
by the following expression, 



M n k n' = J PdS( Z nk\ n ' A*n>% 
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Here e^ost is the host semiconductor dielectric constant, 
\z t) is a six-component envelope-function eigenspinor 
of the Hamiltonian Hq, and the Thomas- Fermi screen- 
ing wavevector was approximated by the parabolic band 
expression, qrF = \J 3e 2 p/ (2ehost£oEF), where p is the 
itinerant hole density and Ef is the Fermi energy. 

Recent experiments have established that magnetic 
and transport properties of (III,Mn)V DMS ferromag- 
nets are sensitive to post-growth annealing protocols 
P,P| Jl^Jl^] , and that this sensitivity is associated with 
changes in the density of defects that compensate the 
Mn 2+ acceptors. Our model recognizes that the trans- 
port properties of these materials are not determined 
solely by the scattering from substitutional Mn 2+ impu- 
rities and allows explicitly for scattering from compen- 
sating defects. We assume that compensation can occur 
due to the presence of As-antisite defects (common in 
low-temperature MBE grown GaAs hosts) or due to Mn 
interstitials. As-antisite defects are non-magnetic and 
contribute only Z = 2 Coulomb scattering. Mn intersti- 
tials, when they are present [|l8|jl6| , are unlikely to be 
magnetically ordered and can also be modeled as Z = 2 
donors p9| . Overall charge neutrality implies that the 
density of holes is p — N Mn 2+ — 2N C , where N c is the 
density of compensating impurities. 

Rough estimates for the T = quasiparticlc scat- 
tering rates can be obtained from parabolic-band- 
approximation expressions for majority heavy-hole states 
assuming Mn 2+ kinetic-exchange scattering only, T p d — 
(iV M „2 + )J 2 d S 2 mV2m*£i7(47r?i 4 ), or Mn 2+ and As- 
antisite Coulomb scattering only which leads to scatter- 
ing rate Tc given by the Brooks-Herring formula [po[ . 
(We show below that most of the current in a (Ga,Mn)As 
ferromagnet is carried by majority spin heavy-holes.) 
Taking the heavy-hole effective mass m* = 0.5m e , p — 
0.4 nm -3 and Mn fraction x — 5%, we obtain UT p d ~ 
20 meV and hTc ~ 150 meV. Our full numerical six- 
band calculations are consistent with these estimates, 
and in particular predict that the Coulomb contribu- 
tion to the elastic scattering rate is several times larger 
than the kinetic- exchange contribution for typical chem- 
ical compositions. The total scattering rate, averaged 



over the majority heavy-hole Fermi surface, decreases 
with increasing density of itinerant holes for a fixed Mn 2+ 
concentration. One important conclusion of this analy- 
sis is that even in the heavily doped and compensated 
(Ga,Mn)As DMSs, the lifetime broadening of the quasi- 
particle (hT) is smaller than the valence band spin-orbit 
coupling strength (A so = 341 meV), which partially jus- 
tifies the neglect of disorder |ll| in evaluating some prop- 
erties, e.g. the anomalous Hall conductivity. 

In Fig.^Jwe plot a xx , calculated numerically using the 
six-band Kohn-Luttinger model and Eqs. (|l|) and (ph, for 
a fully strained Mno.06Gao.94As sample grown on a GaAs 
substrate and assuming compensation due to As-antisites 
alone (i.e. no Mn-interstitials present). The substrate 
- DMS lattice mismatch, eo = (a su b — apMs)/ a DMS, 
is between -0.002 and -0.003 in this case [§Jl6j. (Note 
that a su b and cldms arc the lattice constants of a fully 
relaxed substrate and ferromagnetic layer, respectively.) 
The parameters of the six-band Kohn-Luttinger model 
and strain coefficients used in these calculations are given 
in pi] . We also show in Fig. |] separate contributions from 
individual heavy- and light-hole bands and demonstrate 
that in the ferromagnetic state the current is carried 
mostly by the majority-spin heavy-holes, a property that 
will be important for understanding the spin-injection 
properties of (III,Mn)V DMS ferromagnets. The abso- 
lute conductivities predicted by our model are reason- 
ably consistent with experiment BSJjll. For lower Mn 
concentrations (x < 4%) the theoretical conductivities 
become several times larger than the measured values 
due, we believe, to some combination of inaccuracy in 
our scattering amplitude estimates, sources of disorder 
in current experimental samples that we do not account 
for in the model, and especially at small x coherent scat- 
tering effects that eventually lead to localization. 

As mentioned above, strong spin-orbit coupling in 
the semiconductor valence band leads to a variety 
of magneto-anisotropy effects [Q : ||j2^]. For dc trans- 
port the in-plane conductivity along, e.g. the x- 
direction should change when the magnetization is 
rotated by applying a magnetic field stronger than 
the sample's magneto-crystalline anisotropy field. In 
Fig. H we show the anisotropic magnetoresistance co- 
efficients, AMR lp = [p xx (fl\\x) - p xx (Cl\\y)]/p xx (fl\\y) 



and AAIR 



[p xx (il\\x) - Pxx(tt\\z)]/p xx (fl\\z), for or- 



thogonal magnetization directions in the plane of the 
thin ferromagnetic layer, and for one of the magne- 
tization directions along the growth direction, respec- 
tively. The Mn fraction assuming As-antisite compensa- 
tion alone is indicated in the figure by x\. The total Mn 
fraction (including substitutional Mn 2+ atoms and Mn- 
interstitial atoms) for compensation due to interstitial 
Mn alone is labeled by x%. The plots demonstrate that 
in (Ga,Mn)As/GaAs ferromagnets, which have compres- 
sive strain (eo < 0), AMR is negative for typical chemi- 
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cal compositions, and \AMR ip \ < \AMR op \. (Note that [11] 
AM Ri P = AMR op in unstrained cubic DMSs.) The main 
plot shows that for a fixed hole density (p = 0.4 inn" 3 ) [12] 
the magnitude of the AMR decreases from ~ 10% to P-^l 
~ 1% with increasing Mn fraction. All those observa- 
tions are consistent with available experimental data on 
(Ga,Mn)As DMS's A more detailed compari- 

son to measurements by Gallagher et al. pi shows that 
theoretical data assuming As-antisite compensation only 
overestimate the decrease of AMR with increasing Mn 
fraction. Better quantitative agreement is obtained for 
compensation from Mn-interstitials, which is consistent 
with the presumed dominance of Mn-interstitial defects 
over As-antisite defects in the samples measured by Gal- 
lagher et al. fl6f| . 

In addition to this AMR engineering through dop- 
ing that has been confirmed experimentally our theory 
predicts a large sensitivity of the spontaneous transport 
anisotropy to strain. While strain does not play a signifi- 
cant role for AMRi p , as one might expect from symmetry 
considerations, AMR op can change by more than 10% 
over the range of strains that can be achieved in the thin 
ferromagnetic layers by a proper choice of the substrate 
JjJ. For larger Mn concentrations, AMR op is predicted to 
become positive in samples with tensile strain, as shown 
in the main plot of Fig. ||. 
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FIG. 1. Conductivity for measuring current and magneti- 
zation directed along the x-axis in the plane of the (Ga,Mn) As 
film as a function of the total hole density. These results were 
obtained for a GaAs semiconductor host doped with 6% Mn 
and with strain eo = —0.002. No Mn-interstitial atoms are 
assumed to be present in the ferromagnetic layer. For typical 
hole densities the current is carried mostly by the majority 
heavy-holes. 
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FIG. 2. Anisotropic magnetoresistance coefficients as a 
function of strain (main plot and lower inset) and hole density 
(upper inset). Mn fractions corresponding to compensation 
due to As-antisites alone and Mn-interstitials alone are la- 
beled as x\ and X2 respectively. 
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